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ABSTRACT 

We use deep images taken with the Advanced Camera for Surveys on board the Hubble 
Space Telescope of the disk galaxy NGC 891, to search for globular cluster candidates. 
This galaxy has long been considered to be a close analog in size and structure to the 
Milky Way and is nearly edge-on, facilitating studies of its halo population. These 
extraplanar ACS images, originally intended to study the halo field-star populations, 
reach deep enough to reveal even the faintest globular clusters that would be similar 
to those in the Milky Way. From the three pointings we have identified a total of 
43 candidates after culling by object morphology, magnitude, and colour. We present 
(V, I) photometry for all of these, along with measurements of their effective radius 
and ellipticity. The 16 highest-rank candidates within the whole sample are found to 
fall in very much the same regions of parameter space occupied by the classic Milky 
Way globular clusters. Our provisional conclusion from this survey is that the total 
globular cluster population in NGC 891 as a whole may be almost as large as that of 
the Milky Way. 
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1 INTRODUCTION 

The stellar content, both the diffuse component and star 
clusters, of the outskirts of galaxies are among the old- 
est and the most metal-poor stellar components of galax- 
ies. Their properties are clues to the understanding of how 
galaxies have assembled their mass, and constrain the early 
phases of galaxy formation. Recent Hubble Space Tele- 
scope (HST) imaging and ground-based spectroscopy of star 
clusters have revolutionized our understanding of galaxy 
formation and evolution. However, elliptical and lenticu- 
lar galaxies have received by far the largest amount of 
attention due to their much richer globular cluster (GC) 
systems and freedom from internal extinction. Conversely, 
our knowledge of GC systems in spirals is still limited 
to essentially the Galaxy and M31, supplemented by a 
handful of more distant galaxies ([K isslcr-Pa tig et al. 111999; 
Goudfrooii et al. 1120071; IChandar et al. 1120041 ; iRhode et al. I 
20071 ; ISpitler et al. II2OO6I ; iMora et al. 112007ft . Until very re- 
cently, available data indicated that the GC systems of the 
Galaxy and M31 were quite similar and thus it was nat- 



1 This work was based on observations with the NASA/ESA 
Hubble Space Telescope, obtained at the Space Telescope Science 
Institute, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS 5-26555. 



ural to assume this held true for spirals in general. How- 
ever, in the last few years, a growing body of evidence sug- 
gests that there are some important differences between the 
GCSs of these two key spirals. Among these is the likeli- 
hood that M31 posses ses young, thin disk massive clusters 
|Morrison et al. | 2004), as well as intermediate-age massive 
clusters (e.g. iBeaslev et al. 112004] ). populations that are not 
present amongst the classically old GCs in the Galaxy. The 
young and intermediate-age massive clusters in M31 are 
significantly more massive than any open clusters in the 
Milky Way. Similar clusters ha ve now been seen in a va- 
riety of other spirals (e.g. M33; IChandar et al. 112006ft . The 
GC population in M31 also contains a subset of object of 
extended and di ffuse nature, unlike any c lusters found in 
the Milky Way (|Huxor et al. 1 120051 . 12008ft . These objects 
are found to have simil ar stellar population s to those of the 
Milky Way's old GCs l|Mackev et al. 112006ft , and to fill the 
gap in struc tural-parameter sp ace between GCs and dwarf 
spheroidals l|Huxor et al. 112005ft . 

NGC 891, a nearby large late- type disk galaxy, has often 
been described in shorthand as a "clone" of the Milky Way, 
since it has a very similar total lumin osity, bulge size, an d 
disk with prominent dust lanes (e.g. Ivan der Kruit 1984). 
A stu dy of surface-brightness photometry (|Morrison et al. I 
Il997ft revealed the presence of an extensive thick disk. Its ori- 
entation almost precisely edge-on to our line of sight makes it 
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particularly attractive for studies of its disk and halo stellar 
populations, facilitating comparative studies of components 
such as the thick disk, stellar streams and substructures, 
and the total visible mass and extent of the halo. In addition, 
NGC 891 is close enough to us that HST imaging is easily ca- 
pable of resolving the halo stars, enabling direct star-by-star 
statistical studies of its old stellar populations. Using deep 
HST imaging of three extra-planar fields extending outward 
to more than 10 kpc from the plane of the galaxy, we have 
been studying the resolved stellar populations to investi- 
gate the spatial strucuture, the bi-dimensional distribution, 
to search for substructure, and to constrain the metallicity 
distribution functions. 

One obvious component of its halo that has not been in- 
vestigated to similar detail is the GC system. If this galaxy 
is indeed similar to the Milky Way, then ~ 100 - 200 of 
these luminous, old star clusters should be present and rel- 
atively easy to find. A signal pointing in the opposite direc- 
tion, however, is found in the pion eering effort to search for 
a GC population in NGC 891 by Ivan den Bergh fc Harrisl 
(1982), from starcounts on wide-field photographic plates. 
They found no conclusive evidence for any GCs and an up- 
per limit of Sn < 0.2 on the specific frequency. It should 
be realized, however, that this imaging material was taken 
with l" seeing quality and had a much brighter limiting 
magnitude than is conventionally possible with modern cam- 
eras, making it difficult to find traces of a GC population 
in the presence of significant field contamination. Also, with 
such material, individual GCs cannot be distinguished either 
from foreground stars or faint, small background galaxies, 
and both types of contaminants are present in large num- 
bers within the NGC 891 field. 

In the present paper, for the first time we conduct a 
search for, and characterization of, individual GCs in NGC 
891. Contrary to our original expectations, we find that this 
galaxy does indeed have a roughly normal GC population 
for its size and type. Throughout this paper, we adopt the 
dis tance modulus (m -_M)o = 29.94 {d = 9.7 Mpc) derived 
by iMouhcine et al. 1 l|2007h . along with a foreground red- 
dening E v -i = 0.08 and A v = 0.20 from the NASA/IPAC 
Extragalactic Database (NED). 



2 THE DATABASE AND INITIAL SEARCHES 

The raw data consist of the same deep HST images taken 
with the Wide Field Channel of the Advanced Camera 
for Surveys (A CS/WFC) camera used in previous papers 



in this series ( Mouhcine et al~l 120071 ; llbata et al~l 120081: 
iReikuba et al. I 20081 ) . In these papers, Mouhcine et al. 



2007) discuss the stellar halo population; llbata et al. 



2008) present a structural analysis to establish the pres- 



ence of a thick disk and a stel lar halo as well a s sma ll-scale 
substructures in the halo; and Reikuba et al. I (|2008h inves- 
tigate the stellar populations of the thick disk and the inner 
halo. 

Images were taken at three pointings running parallel 
to, and east of, the disk of NGC 891. A finder chart for 
thes e three fields (labelled HI, H2, H3) is shown as Figure 
1 of IMouhcine et al~l (|2007T l and is reproduced here in Fig- 
ure [1] The three fields were placed to probe both the disk 
and thick disk of the galaxy over a wide radial range as well 



as the inner halo, with slight overlaps between fields to en- 
sure photometric consistency. For each field, the ACS/WFC 
observations consist of three full-orbit integrations in both 
F606W and F814W (for convenience we refer to these below 
as V and I). The photometry reaches to I ~ 29, approxi- 
mately 3 magnitudes below the tip of the red giant branch. 
Artificial-star experiments were performed to estimate the 
completeness of our data set. These were carried out in the 
usual fashion with artificial stars generated from the stellar 
point spread function (PSF) constructed during the pho- 
tometry measurements. The 50 per cent completeness limit, 
for colours typical of old simple stellar populations (i.e., 
(V — /) < 1.3), occurs at I ~ 27. 3 — 28.2 depending on g alac- 
tocentric distance and crowding. IReikuba et al. I (|2008T ) give 
a thorough description of the data, the data reduction pro- 
cess, and the completeness simulations; we therefore refer 
the reader to this paper for the details. 

To start the process of identifying candidate GCs, we 
visually inspected every part of each field and marked any 
objects that might possibly be globular clusters. The selec- 
tion was done independently in each filter. The western sides 
of fields HI and H2 particularly are heavily contaminated 
by the bulge light and planar dust lanes in the galaxy, so to 
aid the identification process we first subtracted off median- 
smoothed versions of the fields (unsharp masking) using a 
39px— square median filter box. The images very obviously 
resolve the stellar halo red-giant population of NGC 891, 
with the brightest normal red giants lying at / ~ 26. The 
faintest known GCs in the Milky Way lie at this level or 
brighter (with the vast majority at Mi < —6, translating 
to I ~ 24 at the distance of NGC 891) and so these images 
should sample virtually the entire GC luminosity distribu- 
tion. The fact that the NGC 891 halo stars are well resolved 
also means that at least some of its GCs could show some 
degree of resolution of their individual red giants outside 
their cores; some examples of these are shown later. 

The main search criteria for this initial round of se- 
lection were that GC candidates should be (a) at least as 
bright as the brightest halo red giants, (b) relatively un- 
crowded, and (c) morphologically symmetric; that is, with 
no features such as small spiral arms, tidal tails, compan- 
ions etc. that would be clear markers of faint background 
galaxies. This initial selection was deliberately made very 
generous; at this stage, no candidate was rejected by ellip- 
ticity, scale size, location, or colour. 



3 REFINING THE CANDIDATE LIST 

With these preliminary lists identified, contaminants were 
weeded out through a series of objective criteria. First was 
to match the lists in both filters and keep only those iden- 
tified on both V and I. This step quickly removes any ex- 
tremely red or extremely blue objects, or potential artifacts 
appearing on one filter only. 

The next stage was to measure the characteristic size 
and shape of each object. W e used the ISHAPE profile- 
fitting code of lLarsen~l jl999) to derive the effective (half- 
light) radius m and ellipticity e = (1 — b/a) of each ob- 
ject. On each image, 15 to 20 moderately bright, uncrowded 
stars were combined with iraf/daophot to construct an em- 
pirical PSF for the frame. For each candidate GC, ISHAPE 
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Figure 1. Location of the three HST/ACS t arget fields along 
the no rtheast side of the disk of NGC 891. See iMouhcine et al, I 
d2007T ) for additional description. 



was then use d to convolve the PSF with a "King30" model 
profile, i.e. a iKing I (|l962l ) model with concentration index 
r t /r c — 30 characteristic of the average for known globu- 
lar clusters. Here as usual r c and rt are the King-model 
core and tidal radii. The assumed effective radius rn and 
axial ratio b/a o f the model a re then varied till a best fit 
is achieved (see lLarsen 1 1 19991 . for exte nsive di s cussio n of 
the technique ). Extensive simulations bv lLarsen~l (Il999l ) and 
lHarris I (|2009l ) show that the derived (rh,e) are nearly in- 
dependent of the assumed r t /r c ratio in situations like this 
one where rn < fwhm(PSF) (see below). 

Because NGC 891 is relatively nearby, the size range 
Th — 2 — 5 parsecs typifying the majority of GCs in the Milky 
Way converts to a range of angular size — 0.04" — 0.11" 
or about 1 to 2.5 pixels on the ACS/WFC camera. Objects 
this extended can be termed "partially resolved" because 
their intrinsic radii are comparable to, or smaller than, the 
PSF, but they are easily distinguishable from stars: The 
stellar PSF on the ACS/WFC has a FWHM = 1.9 ± 0.1 
pixe ls, and extensive tests of the ISHAPE profile fitting code 
(e.g. lLarsen II 19991 ; lHarris 1120091) show that the effective radii 
of partially resolved objects such as these can be correctly 
detected and measured down to 20% of the PSF FWHM, 
and even smaller under conditions of high S/N. Thus all GCs 
with characteristic sizes comparable to those in the Milky 
Way should be easily found. These size measurements can 
then quickly be used to eliminate all stars from our sample 
of candidate GCs. 

Use of the ellipticities is also effective. Small differences 
in mean GC shapes have been found fr om galaxy to galaxy 
that are not yet we ll understood (e.g. lHarris et al. 2002: 
lHan fc Rvden |[l994T ), but these differences are minor, and 
the vast majority of known GCs in any galaxies surveyed 
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Figure 2. Distributions of globular cluster candidates by half- 
light radius and eccentricity, as measured by ISHAPE. Objects 
in NGC 891 are shown as crosses, while the Milky Way globular 
clusters projected to the same distance are shown as solid red 
dots. We reject any candidates lying to the left of or above the 
dashed lines, as being either star like or too elliptical. 



so far are quite round in projected shape (e < 0.15) and 
virtually none are known with e > 0.3. 

In Figure [2] we show the distribution of the candi- 
date objects measured in both V and I by their ISHAPE- 
determined values of and e . For compar ison, the same 
data for the Milky Way GCs (|Harris Ill996l ) are shown as 
well (solid dots in the Figure), where their values are 
projected to the size they would appear if placed at the 9.7- 
Mpc distance of NGC 891. The candidate list has a large 
number of objects at < 0.02"; these are probable stars 
and can rather safely be rejected. We also conservatively 
eliminate any with e > 0.35, a limit well above the most el- 
liptical known Milky Way clusters. Additional visual inspec- 
tion of these highly elongated objects confirms that they are 
likely to be background galaxies; none have any indication 
of being resolved into stars, and many are next to other ob- 
vious galaxies on the fields. The one exception we made to 
this exclusion was for a few very faint candidates that may 
be resolved into stars and were also somewhat elongated. 
Three examples of these are shown in one of the figures be- 
low. These have luminosities and effective radii that, if they 
are indeed clusters, would make them roughly comparable to 
the Palomar-type clusters in the Milky Way or even fainter. 
We regard their true nature as more uncertain than most of 
the rest of the sample, but keep them in the list. 

The third stage of culling is to use colour and magni- 
tude. The raw colour-magnitude diagram (CMD) for the 81 
candidates surviving the structural parameter tests is shown 
in Figure |3] Among these are five which appear on the small 
overlap areas between H1/H2 and H2/H3 and were iden- 
tified on both. In the CMD, we plot the apparent magni- 
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Figure 3. Distribution of candidate GCs in the colour-magnitude 
plane. The NGC 891 objects shown here (open symbols with 
crosses) are ones with effective radii and ellipticities within the 
range of known globular clusters. The Milky Way clusters (solid 
red dots), projected to the distance of NGC 891 and with its fore- 
ground reddening applied, are shown for comparison. The two 
dashed lines show reddening trajectories, with slope A§Q§yy = 
3-E606VK-814W ■ 



tudes and colour indices in the filter system native to ACS 
(F606W, FS14W) , as measured through iraf/daophot and 
aperture photometry with an aperture radius r ap = 3 px 
= 0.15". This aperture size corresponds to about 2.5r;, for a 
median GC in the Milky Way and thus safely includes most 
of its light. Cluster-to-cluster differences in mean that any 
fixed-aperture photometry will not include the same fraction 
of their true total magnitude, but at this stage we are more 
strongly interested in the cluster colors, and for the many 
GC candidates projected on the bulge and disk regions of 
NGC 891, the aperture radius needs to be as small as pos- 
sible to avoid large uncertainties from field contamination. 
The calibration of the photometry follows the recent filter 
zeropoints published by STScI, F606W = 26.420- 2. 51og(/) 
and F8UW = 25.536 - 2.51og(/) where / represents the 
measured counts per second. These magnitudes can be con- 
verted to (V, I ) follo wing the empirical transformations of 
ISirianni et al~l (120051 ). namely V = F606W + 0.236(1/ - I) 
and I = F8UW - 0.002(V - I); we also use these in the 
later discussion. 

The known Milky Way globular clusters (shown for 
comparison in the figure) occupy a narrow range in colour, 
and a range in magnitude that is limited to F606W < 24.5 
if projected to the distance of NGC 891. However, we very 
conservatively define our "best", highest-confidence sample 
of GC candidates to be the ones in the region F606W < 
23.7, (F606W - F8UW) < 1.15. There are 16 such objects, 
which by definition have clearly survived all of our rejection 
tests by scale size, ellipticity, and photometry. Although sev- 



eral others are in the same blue color range and in the fainter 
magnitude range 23.7 < F606W < 24.5, close inspection of 
these fainter ones on the images shows that they are proba- 
bly made up of a wide mixture including small bulge or disk 
clusters, faint and small background galaxies, and perhaps 
the occasional classic GC. In addition, if most of these ac- 
tually were real GCs then it would immediately imply that 
the globular cluster luminosity function (GCLF) in NGC 891 
would be very different from that of the Milky Way, much 
more strongly weighted to the faint end. For these reasons 
we do not rank them in the highest-confidence list. 

We cannot immediately reject candidates that are red- 
der than the Milky Way GC sequence because individual 
objects may be reddened by the heavy dust lanes within 
NGC 891 itself. However, clusters sitting behind significant 
amounts of dust extinction would be most likely to fall 
within the range shown by the upper and lower redden- 
ing lines in the figure. We therefore reject objects falling 
clearly below this region, i.e. those with F606W > 24.5 and 
(F606W - F8UW) < 1.25. Again, close inspection of the 
visual appearance of these very faint objects shows that vir- 
tually all of them are consistent with identification as small, 
distant background galaxies. 

In addition, we can reject very red objects that are lo- 
cated spatially well away from the plane of NGC 891. In 
Figure [4] the colours for the remaining candidates are plot- 
ted against their projected distance above the plane of the 
galaxy. The large circled crosses show the 16 "best" ob- 
jects that fall closest to the normal GC sequence in the 
colour-magnitude diagram. The candidates with (F606W — 
F814W) > 1.17 and Z > 3 kpc can be rejected since they are 
far redder than any known star clusters; note that no mod- 
els of old simple stellar populations can account for such red 
colours, and these same objects are unlikely to be heavily 
reddened since there is no detectable trace of heavy dust ex- 
tinction of any kind at such large distances fro m the plane. 
Using the Hi map of lOosterloo et al. 1 l|2007l ) to estimate 
th e E(B — V) reddening following the Galactic calibration 
of iRachford et al~l l|2008l ) indicates that although the cor- 
rection is substantial near the galactic plane, the correction 
declines rapidly away from the plane, so that at the outer 
edge of the ACS survey region, the estimated internal ex- 
tinction as derived from the Hi column density amounts to 
less than 0.001 mag. 

Our final list of GC candidates, after all the culling 
steps described above, consists of the 43 objects listed in 
Table 1. Successive columns give (1) a running ID number, 
(2) the ACS field on which it lies, (3,4) right ascension and 
declination (J2000) in decimal degree format as calculated 
from the astrometric parameters directly from the image 
headers, (5,6) location (x,y) (pixels) on the particular ACS 
field, (7-8) magnitudes, colours, and internal uncertainties, 
(9) half-light radius rn in arcseconds, (10) ellipse-fitted axial 
ratio (b/a) as measured from ISHAPE, and (11) any com- 
ment on the object ranking or type; here, "LSB" means low 
surface brightness. 

A definitive test of membership in NGC 891 would be 
direct measurement of radial velocity. The mean velocity 
of the galaxy is v r — 528 km s~ , thus its GCs should all 
have 300kms _1 < v r < 800 km s . Because of our morpho- 
logical selection criteria, the final candidate sample has no 
foreground stars, so the only remaining contaminants can be 
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I i i i i I i i i i I i H Figure 5. Six of the candidate GCs in NGC 891, showing their 

5 10 marginal resolution into stars. From Table 1, these are G13, G24, 

Z (kpc) and G34 (top row) and G22, G15, and G28 (bottom row). Each 

image here is 100 px across (5"). Note the large numbers of field 
halo stars within NGC 891 across each field. 



Figure 4. Measured colour of candidates versus their distance 
above the plane of NGC 891. The 16 objects shown as large circled 
crosses make up the "best" candidate GC sample that are closest 
to the normal Milky Way cluster sequence in the previous figure. 
All other, lower-quality candidates are plotted as crosses. The 
objects to the upper right of the dashed lines are red ones at large 
distances from the plane of the galaxy and can be rejected (see 
text). Note that two objects just on the borderlines are included 
in the final candidate list. 



faint background galaxies and perhaps a few massive open 
clusters in the disk or bulge of NGC 891 itself. However, 
only 14 candidates are brighter than V ~ 23, making veloc- 
ity measurements challenging for all but these few. 



4 DISCUSSION 

After the many stages of weeding out individual contami- 
nants, it is encouraging that we have numerous objects that 
do indeed resemble normal GCs closely even though the var- 
ious culling stages were quite conservative. In Figure [5] we 
show thumbnail images of six of the candidates, demonstrat- 
ing their partial resolution into stars. In the top row, three of 
the brightest are shown, and in the bottom row, three of the 
faintest. As noted above, we regard the identifications of the 
faintest ones as generally more uncertain, and in some cases 
it may be more likely that they are actually background 
galaxies. However, we prefer to keep some contaminants in 
the list rather than to miss a few real clusters. 

More quantitative comparisons can be seen in Figures 
E] and [3 These two fi gures illustrate that our 16 best candi- 
dates have colours, magnitudes, half-light radii, and elliptic- 
ities all closely resembling those of normal GCs in the Milky 
Way. Full histogram comparisons between the two galax- 
ies to search for any finer differences would be inconclusive 
given the small sample size of NGC 891 GCs. The lack of 
NGC 891 candidates with formal ellipticities smaller than 
~ 0.05 (compared with the Milky Way, which has many 
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Figure 6. The colour-magnitude diagram for the final sample of 
43 GC candidates, compared with the Milky Way clusters (solid 
red dots) as before. The 16 best candidates are plotted as cir- 
cled crosses and the remaining 27 lower-probability candidates as 
smaller crosses. 



such clusters) cannot be assigned too much weight given 
the internal fitting uncertainties in ISHAPE for these par- 
tially resolved objects; but if real, it would indicate that 
the NGC 891 clusters are more similar i n ellipticity to those 
in the Magellanic Cl ouds or NGC 5128 (|Harris et al. 1120021 ; 
lHan fc Rvd"en~lll994r i. 

The plot of th versus luminosity shows hints of the same 
trends in the Milky Way clusters, namely that the lower 
envelope of points shows gradually increasing cluster scale 
size with luminosity; and that the scatter in ru increases at 
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Figure 7. Ellipticity (upper panel) and half-light radius (lower 
panel) plotted versus cluster luminosity. As before, the Milky Way 
clusters are plotted as solid red dots, the 16 best NGC 891 can- 
didates as large circled crosses, and the 27 other candidates as 
smaller crosses. 



Figure 9. Cluster effective radius plotted versus projected 
galactocentric distance Rgc- The large circled crosses show the 
16 best candidates from Table 1 while the smaller crosses are for 
the remaining less certain candidates. As in previous diagrams, 
the red solid dots show the same data for the Milky Way clusters. 
Most of the largest candidates are also the faintest, thus their 
identification as true clusters is more uncertain; see text. 
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Figure 8. Spatial locations of the selected GC candidates, super- 
imposed on the positions of the Milky Way clusters to the same 
scale. A schematic outline of the Milky Way bulge and disk are 
drawn in for scale comparison. In this plane, the A"— axis points 
roughly northward and Z westward. The 16 best GC candidates 
in NGC 891 are plotted as large red dots, the remaining 27 as 
small blue dots, and the reference Milky Way clusters as small 
crosses. 



fainter luminosity. The GC candidate with the biggest scale 
size, G40 at rn — 24 pc (lying at the uppermost edge of 
the lower panel in Fig. [7}, is indeed an extended and bright 
object but does not show strong traces of resolution into 
stars and may be a background galaxy. Deeper images or 
radial velocity measurement would make a definitive test. 

The other 27 candidates mostly occupy the low- 
luminosity sides of both graphs, and show larger dispersions 
in both e and r h that also match the Milky Way spread quite 
well. Unfortunately, however, many of these candidates are 
likely to be background galaxies as well. 

As another comparison with the Milky Way, we show 
the spatial locations of the NGC 891 clusters in Figure [5] 
where (Y, Z) in kiloparsec units are the projected distance 
components parallel and perpendicular to the plane of the 
galaxy. The galactic center is at (0,0). The orientation o f 
the figure is chosen to match Figure 3 in llbata et al. I (2008). 
Here, as a schematic comparison, the positions of the Milky 
Way GCs are shown projected on the same plane, essen- 
tially viewing our galaxy as it would be seen from outside, 
and looking inward along the line connecting the Sun to the 
Galactic center. This choice of two-dimensional projection 
minimizes the effects of random distance errors to individual 
Milky Way clusters, most of which lie toward the Galactic 
center. 

Lastly, in Figure [9] the distribution of cluster scale size 
rh versus galactocentric distance is shown. In this case the 
Milky Way data (small dots) are for projected Galactocen- 
tric distance Rgc = VY 2 + Z 2 to make the graph for the 
two galaxies strictly comparable. In the Milky Way we see 
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the well known tre nd for the mean r^ to increase gradu- 
ally w ith Rgc (e.g. Ivan den Bergh aL~lll99ll ; [Jordan et al. I 
120051 ) and for the range in rh to increase at larger Rgc as 
well. The NGC 891 data hint barely at the same increase 
amidst the large scatter. Five of the NGC 891 clusters have 
radii that fall clearly above the Milky Way distribution, how- 
ever; these are G07, G08, G13, G19, and G40. If at least 
some of these five are not contaminants, we can speculate 
that they may have originated in dwarf satellite galaxies 
within which GCs ha ve characteristically larger radii (see 
|Pa Costa et al. II2009I . for a recent discussion). Such an ori- 
gin w ould mesh well with our analysis of the NGC 891 stellar 
halo l|lbata et al. I I2008T ) which shows significant substruc- 
ture and indicates that at least some part of it has been 
accreted from satellites. Without more material than our 
presently rather slim GC sample to base it on, it is risky to 
carry such discussion further. 

Field HI is closest to the galaxy center and bears out 
our expectations that it should contain the most GCs of the 
three fields (see Fig. [Sj . Furthermore, any GCs particularly 
in HI that happen to lie behind the heavy dust lanes there 
may well be escaping detection, so by symmetry we might 
then expect that the true total GC population in HI is as 
much as twice as large as the 10 good candidates and 9 
lesser-quality ones we actually see. Combining all the pre- 
ceding arguments, we estimate very roughly that our three 
observed fields may contain 25 ± 5 high-quality GC candi- 
dates over all magnitudes after correction for contamination 
and for losses due to extinction. 

Estimating the total GC population within the entire 
galaxy requires a large extrapolation from the severely lim- 
ited area coverage of our survey. Our three ACS fields cover 
the region to the east side of the northern disk, so they 
can take in at most ~ 1/4 of the entire population. We 
therefore suggest, albeit tentatively, that the total should 
be Nt > 100 and perhaps as large as ~ 140. Alternately, 
we can compare the number of cluster candidates in our 
three fields with the number of Milky Way clusters that fall 
within the same area when projected onto the same scale 
(see Fig. |8j). In the Figure, 30 to 35 Milky Way clusters fall 
within the marked-out areas of Hl/2/3, compared with any- 
where from 16 to 43 candidates in our defined list for NGC 
891. Since the total cataloged population of Milky Way clus- 
ters is approximately 150, these ratios suggest that the total 
population of NGC 891 globular clusters is in the range of 
~ 80 to 200, in good agreement w ith the previous estimate. 
The resulting specific frequency (Ivan den Bergh fc Harris! 
1 19821 ) is then Sn ~ 0.3 — 0.5, if we adopt a galaxy luminosity 
Vt = 8.82 from the NED database. This specific frequency 
would put NGC 891 reasonably in the range of most know n 
disk galaxies similar to the Milky Way (e.g. lHarris ll200ll ). 

In summary, our admittedly restricted look at the glob- 
ular cluster system of NGC 891 leads us to conclude that 
it resembles that of the Milky Way rather closely in ap- 
proximate total numbers, spatial distribution, and struc- 
tural properties of the individual clusters. Our data clearly 
cover only a very limited region of the NGC 891 halo, and 
a more comprehensive targetted search should be able to 
increase the GC sample by a factor of five or more. 
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Table 1. Candidate Globular Clusters in NGC 891 



ID Field RA Dec x (px) y (px) F606W (F606W - F8UW) r h (arcscc) (b/a) Comment 



G01 


HI 


35 


.624222 


42. 


.338615 


5480. 


.3 


3664.0 


22 


.089 


G02 


HI 


35 


.625717 


42 


.338867 


5400. 


.6 


3682.1 


24 


.065 


G03 


HI 


35 


.625805 


42. 


339149 


5396 





3702.4 


25 


.020 


G04 


HI 


35 


.627125 


42. 


.339451 


5325. 


,7 


3724.1 


24 


.245 


G05 


HI 


35 


.627289 


42 


.340481 


5316 


.9 


3798.3 


25 


.323 


G06 


HI 


35 


.629326 


42 


.331100 


5209. 





3122.7 


23 


.243 


G07 


HI 


35 


.629837 


42. 


.335514 


5181. 


,5 


3440.6 


23 


.077 


G08 


HI 


35 


.631474 


42. 


330360 


5094. 


,6 


3069.6 


23 


.366 


G09 


HI 


35 


.631748 


42. 


332382 


5079. 


9 


3215.0 


23 


.669 


G10 


HI 


35 


.634480 


42. 


.338356 


4934. 


3 


3645.2 


23 


.996 


Gil 


HI 


35 


.635334 


42. 


.338711 


4888. 


.8 


3670.6 


22 


.981 


G12 


HI 


35 


.636715 


42. 


.336533 


4815. 


4 


3513.9 


24 


.476 


G13 


HI 


35 


.640335 


42. 


.342419 


4622. 


.6 


3937.4 


21 


.301 


G14 


HI 


35 


.640930 


42. 


.327168 


4591. 


.4 


2839.4 


20 


.320 


G15 


HI 


35 


.647358 


42. 


.309727 


4249. 


,5 


1583.5 


25 


.388 


G16 


HI 


35 


.648632 


42. 


320786 


4181. 


,6 


2379.6 


24 


.399 


G17 


HI 


35 


.650364 


42 


.315956 


4089. 


.4 


2031.9 


22 


.673 


G18 


HI 


35 


.650578 


42. 


.340332 


4077. 


.4 


3786.9 


22 


.450 


G19 


H2 


35 


.668228 


42 


.389709 


4683. 


.5 


3592.4 


23 


.215 


G20 


H2 


35 


.668983 


42. 


.389130 


4643. 


3 


3550.5 


23 


.872 


G21 


H2 


35 


.669926 


42. 


.379181 


4593. 


.4 


2834.3 


23 


.140 


G22 


HI 


35 


.670647 


42. 


.343815 


3009. 


.4 


4037.5 


24 


.646 


G23 


HI 


35 


.672882 


42. 


.338108 


2890. 


.5 


3626.5 


22 


.464 


G24 


H2 


35 


.678722 


42. 


396378 


4125. 


4 


4072.4 


21 


.355 


G25 


H2 


35 


.679241 


42. 


.392036 


4097. 


.8 


3759.7 


23 


.972 


G26 


H2 


35 


.682297 


42 


.399582 


3935. 


3 


4302.9 


23 


.150 


G27 


H2 


35 


.684418 


42. 


.401455 


3822. 


.5 


4437.9 


24 


.043 


G28 


H2 


35 


.684708 


42. 


.366974 


3807. 


,5 


1955.0 


25 


.844 


G29 


H2 


35 


.685135 


42. 


.401295 


3784. 


,2 


4426.2 


23 


.698 


G30 


H2 


35 


.686897 


42. 


.400967 


3690. 


.6 


4402.6 


23 


.132 


G31 


HI 


35 


.688915 


42 


.326588 


2037. 


.0 


2797.1 


24 


.181 


G32 


H3 


35 


.689495 


42. 


.432602 


5072. 


.7 


2963.6 


24 


.260 


G33 


H2 


35 


.689682 


42. 


.348801 


3542. 


9 


646.6 


23 


.897 


G34 


H3 


35 


.692303 


42. 


.428921 


4923. 


.5 


2698.4 


21 


.872 


G35 


HI 


35 


.693783 


42. 


.332630 


1778. 


.0 


3232.2 


24 


.057 


G36 


H3 


35 


.696560 


42 


.433994 


4697. 


.0 


3063.5 


25 


.059 


G37 


H3 


35 


.698162 


42. 


.433067 


4612. 


,1 


2996.9 


24 


.165 


G38 


H3 


35 


.701077 


42. 


.428608 


4457. 


.4 


2675.5 


21 


.846 


G39 


H2 


35 


.701416 


42. 


.400101 


2918. 


.6 


4340.1 


21 


.651 


G40 


H2 


35 


.719254 


42. 


.415409 


1970. 


6 


5442.5 


22 


.888 


G41 


H2 


35 


.727627 


42. 


.381878 


1524. 


,6 


3028.4 


24 


.070 


G42 


H2 


35 


.736752 


42 


.386242 


1039. 


,4 


3342.8 


21 


.844 


G43 


H2 


35 


.740395 


42. 


.379120 


845.3 


2830.1 


24 


,223 



± 


0.007 


0.811 


± 


0.013 


0.0342 


0.91 


best 


± 


0.041 


1.126 


± 


0.070 


0.0828 


0.79 


bulge 


± 


0.086 


2.344 


± 


0.093 


0.0397 


0.83 


bulge 


± 


0.111 


0.522 


± 


0.216 


0.1530 


0.70 


bulge 


± 


0.291 


1.806 


± 


0.319 


0.2844 


0.84 


bulge 


± 


0.035 


1.289 


± 


0.048 


0.1025 


0.82 


disk 


± 


0.031 


0.990 


± 


0.051 


0.1576 


0.85 


best 


± 


0.029 


0.814 


± 


0.050 


0.1939 


0.88 


best 


± 


0.053 


1.054 


± 


0.084 


0.0217 


0.76 


best 


± 


0.035 


1.046 


± 


0.068 


0.0222 


0.83 


disk 


± 


0.012 


0.889 


± 


0.031 


0.0250 


0.77 


best 


± 


0.143 


1.928 


± 


0.155 


0.0246 


0.71 


disk 


± 


0.006 


0.923 


± 


0.009 


0.1609 


0.94 


best 


± 


0.002 


0.871 


± 


0.003 


0.0440 


0.94 


best 


± 


0.034 


1.565 


± 


0.039 


0.1142 


0.88 


LSB 


± 


0.018 


1.258 


± 


0.023 


0.4173 


0.36 


galaxy? 


± 


0.005 


0.900 


± 


0.007 


0.0666 


0.91 


best 


± 


0.006 


0.789 


± 


0.009 


0.0669 


0.81 


best 


± 


0.014 


0.946 


± 


0.021 


0.2499 


0.90 


best 


± 


0.020 


1.083 


± 


0.032 


0.2574 


0.73 


disk 


± 


0.007 


1.581 


± 


0.009 


0.1266 


0.84 


galaxy? 


± 


0.023 


0.823 


± 


0.033 


0.4319 


0.69 


LSB 


± 


0.004 


0.681 


± 


0.006 


0.0321 


0.91 


best 


± 


0.003 


0.909 


± 


0.004 


0.0573 


0.83 


best 


± 


0.014 


1.342 


± 


0.018 


0.3257 


0.43 


galaxy? 


± 


0.006 


2.318 


± 


0.007 


0.0285 


0.87 


cluster? 


± 


0.011 


0.925 


± 


0.015 


0.3288 


0.81 


cluster? 


± 


0.044 


1.008 


± 


0.059 


0.4078 


0.36 


LSB 


± 


0.008 


1.727 


± 


0.009 


0.0202 


0.82 


star? 


± 


0.006 


2.201 


± 


0.007 


0.0238 


0.70 


star? 


± 


0.010 


0.825 


± 


0.015 


0.1472 


0.79 


galaxy? 


± 


0.015 


0.695 


± 


0.024 


0.2386 


0.86 


galaxy? 


± 


0.009 


0.638 


± 


0.013 


0.2577 


0.68 


galaxy? 


± 


0.003 


0.925 


± 


0.004 


0.0680 


0.85 


best 


± 


0.010 


0.672 


± 


0.015 


0.0479 


0.74 


galaxy? 


± 


0.022 


2.172 


± 


0.024 


0.1496 


0.75 


galaxy? 


± 


0.011 


1.165 


± 


0.014 


0.2167 


0.91 


galaxy? 


± 


0.003 


0.821 


± 


0.004 


0.1177 


0.94 


best 


± 


0.003 


0.811 


± 


0.004 


0.0685 


0.91 


best 


± 


0.005 


0.813 


± 


0.007 


0.5212 


0.88 


best 


± 


0.010 


0.636 


± 


0.015 


0.1566 


0.71 


galaxy? 


± 


0.003 


1.190 


± 


0.004 


0.0599 


0.83 


galaxy? 


± 


0.010 


0.631 


± 


0.016 


0.2842 


0.90 


galaxy? 
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